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We compute the matrix elements of the electromagnetic (EM) operator between kaon and pion 
states, using lattice QCD with maximally twisted-mass fermions and two flavors of dynamical 
quarks (Nf — 2). The EM operator is renormalized non-perturbatively in the RI'/MOM scheme 
and our simulations cover pion masses as light as 260 MeV and three values of the lattice spacing, 
ranging from ^ 0.07 up to ^ 0. 1 fm. At the physical point our preliminary result for the K ^ n 
tensor form factor at zero-momentum transfer is fj^{Q) = 0.42(2stat), which differs significantly 
from the old quenched result /f'^(O) = 0.78(6) obtained by the SPQcdR Collaboration [|l| with 
pion masses above 500 MeV. We investigate the source of this difference and conclude that it is 
mainly related to the chiral extrapolation of the quenched data. For the case of the tensor charge 
of the pion we obtain the preliminary value fj^{Q) = 0.200(14stat), which can be compared with 
the result /f''(0) = 0.216(34) obtained at A?/ = 2 by the QCDSF Collaboration [|] using higher 
pion masses. 
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1. Introduction 

Precision measurements of weak decays can constrain the parameters of the Standard Model 
(SM) and place bounds on new physics (NP) models, such as supersymmetry. In particular, penguin 
operators between kaon and pion states can place strong bounds on the CP- violating parameters in 
the light quark sector. 

In this contribution we present a lattice study of the electromagnetic (EM) operator, relevant 
in the CP violating part of the K — ni^i— semileptonic decays, performed using the gauge config- 
urations generated by the European Twisted Mass Collaboration (ETMC) with Nf = 2 maximally 
twisted-mass fermions. The EM operator involved in the weak s ^ d transition is given by 

QEM = sF^^a^''d , (1.1) 

where F^y is the EM field tensor. Therefore its matrix element between kaon and pion states 
involves the one of the weak tensor current, which can be written in terms of a single form factor, 
/r(^'),as 

{n\p,)\sa^^d\K\pK)) = {pU-pIp'k) - (1-2) 

where = {pK — Pn)^ is the 4-momentum transfer. Note that the mass factor {Mk + M-^)^^ is 



conventionally inserted in Eq. ( |1.2[ ) in order to make the tensor form factor dimensionless. 

Our simulations cover pion masses as light as 260 MeV and three values of the lattice spacing, 
ranging from ~ 0.07 up to ~ 0. 1 fm. At the physical point our preliminary result for the K ^ n 
tensor form factor at zero-momentum transfer is 

/f ^(0) = 0.42 (2stat) (ETMC) , (1.3) 



where the error is statistical only. Our finding ( |1.3| ) differs significantly from the old quenched 



result fj^'iO) = 0.78(6) obtained in Ref. ^] by the SPQcdR Collaboration with pion masses above 
~ 500 MeV. The reason is mainly due to the non-analytic behavior of the tensor form factor fj^{0) 
in terms of the quark masses introduced by the mass factor {Mk+Mj^)^^ in the parameterization 
(L2). Such a behavior was not taken into account in Ref. [|l|] (see later on). 

In the case of the degenerate tt — )• tt transition, making use of the pion mass dependence 
predicted by Chiral Perturbation Theory (ChPT) in Ref. [^, we obtain for the tensor form factor 
fx^{0), also known as the tensor charge of the pion, the preliminary value at the physical point 

/?'^(0) = 0.200 (14stat) (ETMC) , (1.4) 

which can be compared with the result fj^{0) = 0.216(34) obtained at Nf = 2 by the QCDSF 
Collaboration using simulations at higher pion masses. 

2. K-^ It results 

We have performed the calculations of all the relevant 2-point and 3-point correlation functions 
using the ETMC gauge configurations with Nf = 2 dynamical twisted-mass quarks [Q] generated 



2 



Electromagnetic operator between kaon and pion states 



I. Baum 



at three values of j3, namely the ensembles A2 — A4 at j8 = 3.8 (a ~ 0.103 fm), Bi — B7 at j8 = 3.9 
{a ~ 0.088 fm), and Ci — C3 at j8 = 4.05 (a ~ 0.070 fm). The pion mass ranges between 
~ 260 MeV and ~ 575 MeV and the size L of our lattices guarantees that Mj^L is larger than 
~ 3.3. For each pion mass and lattice spacing we have used several values of the (bare) strange 
quark mass to allow for a smooth, local interpolation of our results to the physical value of 
rus (see Ref. [^]). The calculation of the 2- and 3-point correlation functions has been carried 
out using all-to-all quark propagators evaluated with the one-end-trick stochastic procedure and 
adopting non-periodic boundary conditions which make arbitrarily small momenta accessible. All 
the necessary formulae can be easily inferred from Ref. [^, where the degenerate case of the vector 
pion form factor is illustrated in details. For each pion mass the statistical errors are evaluated with 
the jackknife procedure. 

The tensor current was renormahzed non-perturbatively in the RI'/MOM scheme in Ref. [Q], 
including 0'{a^) corrections coming from lattice perturbation theory [^J. The numerical values used 
in our analyses for the tensor renormalization constant are Z7-(MS,2 GeV) = 0.733(9), 0.743(5), 
0.777(6) for jS = 3.8, 3.9, 4.05, respectively. 

At each pion and kaon masses we determine the tensor form factor fj^{q^) for several values 
of < ql^^^ = {Mk —MjiY in order to interpolate at q^ = 0. Note that, because of the vanishing 
of the Lorentz structure in Eq. (L2), it is not possible to determine fj^{q^) at q^ = q^nax- this 
respect we take advantage of the non-periodic boundary conditions to reach values of q^ quite close 
to q^ = 0. The momentum dependence of fj^{q^) can be nicely fitted either by a pole behavior 

/f"(^')=/f"(0)/(l-4"^') (2.1) 

or by a quadratic fit in q^ 

/r(^')=/f"(0)-(l+4"'?' + cf 9') . (2.2) 

The good quality of both fits is illustrated in Fig. |l], where the results obtained at two different 
lattice volumes are also compared. It can clearly be seen that: i) finite size effects are well below 
the statistical precision of our lattice points; ii) the results for fj^ (0) (as well as those for the slope 
Sj^), obtained using the pole dominance (2A), differ only very slightly from those obtained via the 
quadratic fit (2.2). Therefore in what follows we limit ourselves to the results for both fj^{0) and 



the slope Sj''^ obtained through the pole fit ( |2.1[ ) 

The values obtained for /^'^ (0) and Sj''^ depend on both pion and kaon masses. The depen- 
dence on the latter is shown in Fig. § for fj^{0) at Mj^ ~ 435 MeV and it appears to be quite 
smooth. Thus an interpolation at the physical strange quark mass can be easily performed using 
quadratic splines. This is obtained by fixing the combination (2M|^ — M|) at its physical value, 
which at each pion mass defines a reference kaon mass, M'^^ , given by 

2[M''^ff -Ml = l\Mf^'f - \Mp'f (2.3) 

with Uf^'" = 135.0 MeV and Mp' = 494.4 MeV. 

The results for fj^{0) and the slope Sj^, interpolated at the reference kaon mass Mk = M^k^ , 
are shown in Figs. ^ and Q respectively, for the three lattice spacings of our simulations. It can 
clearly be seen that discretization effects are sub-dominant and therefore, in what follows, we 
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Figure 1: The tensor form factor fj^{q^) obtained at Myc ~ 300 MeV and Mk ~ 530 MeV versus in 
physical units. The dots and the squares (shifted for better clarity) correspond to the gauge ensembles Bi 
and B-i, respectively, which differs only for the lattice size. The solid and dashed hnes are the results of the 
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Figure 2: Results for fj^ (0) versus M\ at ~ 435 MeV. The square corresponds to the value of fj^ (0) 



obtained by local interpolation via quadratic splines (dotted line) at the reference kaon mass M''^f ~ 575 



MeV from Eq. (2.3) 



concentrate on the chiral extrapolation of our lattice data, which in the case of fj^{0) is a much 
more delicate point. 

In Ref. [|l|] the first lattice calculation of the EM operator matrix element between kaon and 
pion states was carried out in the quenched approximation and for pion masses above ~ 500 MeV. 
There the chiral extrapolation was performed assuming that fj'"' (0) reaches a non- vanishing value 
in the SU(3) chiral limit {Mk, Mj^) — )• 0. A simple linear fit in the squared kaon and pion masses 
was attempted obtaining at the physical point the result f^^{0) = 0.78(6). 

In the degenerate case = the chiral expansion of the tensor current has been studied in 
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Figure 3: Results for /f''(0) versus at Mk = M'^-^ in physical units. The dots, squares and triangles 
are our results for the three lattice spacings of the ETMC simulations, specified in the inset. The solid and 



dashed lines correspond to the fit given by Eq. (2.4) with B = and D — Q, respectively. 
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Figure 4: The same as in Fig. ||, but for the slope Sj^ of the tensor form factor at ~ 0. 



Ref. ig]. The main finding is that the form factor fj^ (0) vanishes like Mji for M,^ — )• 0, so that the 
ratio fj^{0)/Mjt tends to a non- vanishing value in the chiral limit. The same argument is expected 
to hold as well in the case of the K ^ n transition: the form factor ff^{0) must vanish in the SU(3) 
chiral Umit in such a way that the ratio ff^{0) / {Mk+Mji) does not vanish in limit Mk = Mji = 0. 
Therefore we perform the chiral extrapolation of our lattice data using the ansatz 



/f '^(O) = (m7 +M^) a [1 +BMl \og{Ml) +CMl + DM^^ , 



(2.4) 



where A, B, C and D are unknown low-energy constants (LECs). The results of the fit (2.4) as- 
suming either B = (no chiral logs) or D = are shown in Fig. ^ by the solid and dashed lines, 
respectively. It can be seen that the effects of the chiral logs are not visible in our data. At the 
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physical point we get 

/f'^(0)=0.42(2s,at) (ETMC), (2.5) 



where the error is statistical only. Had we neglected the mass factor {M'^^ +Mn) in Eq. the 
result at the physical point would change only marginally: fT^{0) = 0.46(2stat)- 

On the contrary, in the case of the quenched data of Ref. [|l]], which were determined at pion 
masses above ~ 500 MeV, the inclusion of the mass factor {Mk+Mj^) in the chiral extrapolation 
changes significantly the result at the physical point by many standard deviations, namely from 
/|''^(0)= 0.78(6) to 

/f^(0) =0.49(4) (SPQcdR). (2.6) 

These findings indicate that the real effect of the quenched approximation does not exceed 15%, 
provided the correct mass factor is included in the chiral extrapolation. 

In the case of the slope Sj'"^ no mass factor should be considered and the chiral extrapolation of 
the lattice data shown in Fig. ^ provides at the physical point the value Sj^ = 1.29(18stat) GeV^^, 
which is consistent within the errors with the quenched result of Ref. ^s^"" = l.n{5) GeW^. 

3. Pion tensor charge 

Following Ref. [||] the chiral expansion of the pion tensor charge fj^ (0) has the form 



M: 



2 



1 + (4^ +C'Ml + D'Ml 



(3.1) 



where /;t ~ 130 MeV and the presence of the mass factor Mji is expected to have an important, 
bending effect on the value extrapolated at the physical point. 

Our results for /^''(O), obtained at jS = 3.9 for 260 MeV < M;^ < 575 MeV, are shown in 
Fig. H and compared with the ones from Ref. [||], having Mji > 440 MeV, and with the quenched 
calculations of Ref. [[T]], ranging from M^^ ~ 530 MeV up to Mj^ ~ 800 MeV. It can be seen that our 
results have a better statistical precision and cover much lighter pion masses, where the bending 
effect due to the overall mass factor Mj^ is clearly visible. 



Using Eq. ( p.l\ ) with our lattice points, we get at the physical point 

f^'^iO) = 0.200 (14,tat) (ETMC) , (3.2) 

to be compared with the QCDSF result [|] at A^/ = 2 

/^'^(O)l = 0.216 (34) (QCDSF) . (3.3) 

Finally, we also apply a simple fit of the form fj^(0) = Mji A' [l +C'MI] to the three quenched 
data of Ref. ||l|], obtaining at the physical point the result 

/^^(O) = 0.221 (21stat) (SPQcdR) , (3.4) 

which clearly shows that quenching effects are sub-dominant on the pion tensor charge. 
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Figure 5: Results for the pion tensor charge fj^ifi) versus in physical units from our simulations at 
j3 — 3.9 (dots) and from Refs. (triangles) and (squares). The solid line is the result of the fit (3.1) 
applied to our lattice points, while the dotted line corresponds to the fit described in the text and applied to 
the quenched data of Ref. 



References 

[1] D. Becirevic, V. Lubicz, G. Martinelli and F. Mescia [SPQcdR Collaboration], Phys. Lett. B 501, 98 
(2001) [ai-Xiv:hep-ph/0010349]. 

[2] D. Brommel et al. [QCDSF and UKQCD Collaborations], Phys. Rev. Lett. 101, 122001 (2008) 
[arXiv:0708.2249 [hep-lat]]. 

[3] M. Diehl, A. Manashov, A. Schafer, Phys. Lett. B622 (2005) 69-82 [hep-ph/0505269]; Eur Phys. J. 
A31 (2007) 335-355 [hep-ph/06 11101]. 

[4] R. Baron et al. [ETM Collaboration], JHEP 1008 (2010) 097 [arXiv:091 1.5061 [hep-lat]]. 

Ph. Boucaud et al. [ETM Collaboration], Phys. Lett. B 650 (2007) 304 [arXiv:hep-lat/0701012]. 

Ph. Boucaud et al. [ETM Collaboration], Comput. Phys. Commun. 179 (2008) 695 [arXiv:0803.0224 
[hep-lat]]. 

[5] B. Blossier et al. [ETM Collaboration], Phys. Rev. D82 (2010) 1 14513 [arXiv: 1010.3659 [hep-lat]]; 
PoS LATTICE2010 (2010)239 [arXiv: 101 1.1 862 [hep-lat]]. 

[6] R. Frezzotti ef aZ. [ETM Collaboration], Phys. Rev. D79 (2009) 074506 [arXiv:0812.4042 [hep-lat]]. 

[7] M. Constantinou et al. [ETM Collaboration], JHEP 1008, 068 (2010) [arXiv: 1004.1 1 15 [hep-lat]]. 

[8] M. Constantinou, V. Lubicz, H. Panagopoulos and F Stylianou, JHEP 0910, 064 (2009) 
[arXiv:0907.0381 [hep-lat]]. 



7 



